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Abstract 

We consider the scaling Lee- Yang model. It corresponds to the unique perturbation 
of the minimal CFT model M2,5. This is not a unitary model. We used known 
expression for form factors in order to obtain a closed expression for a correlation 
function of a trace of energy-momentum tensor. This expression is a determinant of 
an integral operator. Similar determinant representation were proven to be useful 
not only for quantum correlation functions but also in matrix models. 
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1 Introduction 



The theory of massive, relativistic, integrable models is an important part of modern 
quantum field theory Scattering matrices in these models factorize into a 



product of two-body S-matrices |2lf. Form factors can be calculated on the basis of a 



bootstrap approach pl|— 128 



The purpose of this paper is to calculate correlation functions. As usual correlation 
function can be represented as an infinite series of form factors contributions. In this 
paper we sum up all these contributions and obtain a closed expression for the correlation 
function of the energy- momentum tensor. We follow an approach of the paper ||. We 
introduce an auxiliary Fock space and auxiliary Bose fields (we shall call them dual fields). 
These fields help us to represent the form factor decomposition of a correlation function 
in a form similar to "free fermionic" case. This approach was developed in |2^, |10[ . 
Finally a correlation function is represented as a vacuum mean value (in the auxiliary 
Fock space) of a determinant of an integral operator ( |4.3|) . Determinant representation 
were proven to be useful not only for quantum correlation functions || - [[l(J but also 



in matrix models JH]] - [2D]. It helps for asymptotical analysis of quantum correlation 



functions. Among other things this approach helped to calculate asymptotic of time and 



temperature dependent correlation function in Nonlinear Schrodinger equation |L3 . 

The scaling Lee- Yang model can be described by the 01^-perturbation of the non- 
unitary minimal model il^s [[[]]. The theory is known to be integrable and is described 
by factorizable scattering theory of only one kind of particle with mass m. The two-body 
scattering amplitude is given by ]]J 

= sinh/3 + isin(27r/3) 
{P) sinh/3-2sin(27r/3)' 1 ' ' 

The pole at ft — 2ni/3 corresponds to a bound state pole. The residue of the pole is 
negative. So three point coupling is imaginary. These facts come from non-unitarity 
of the theory. This S-matrix can be obtained from breather-breather S-matrix of Sine 
Gordon M. The breather-breather S-matrix of Sine Gordon was first calculated in ||. 

Form factors for the trace of the energy- momentum tensor = T£/4 is defined by 
the matrix element between the vacuum state (vac| and n particle states characterized 
by their rapidities fti (i — 1, . . . , n): 

F n (fti, ...,ftn) = (vac|e(0)|A, ...,ft n ). (1.2) 
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The multi-particle form factors F n was determined in [0, 



F n (f3 u . . . , AO = H n Q n (x u . . . , x n ) I] — , (1-3) 



where Xj = e 13 ^ i = 1 n 



!•••)'"! 



Tim 



where 



H " = -WA^¥)) ■ (L4) 

The function /(/3) is given by 

fw-z&jr^** (i.5) 

The function /(/3) has a single pole in the strip < lm(3 < tt at (5 = 27ri/3 and a single 
zero at (5 = 0. 

We shall use the elementary symmetric polynomials <Jk{xi, ■ ■ ■ , %n), 

n 

Y[(x + Xj) = x n ' k a k (xi, . . .,x n ), 
j=i kez 

It is important to remember that is equal to zero if k < or k > n . The symmetric 
polynomials Q n (xx, . . . , x n ) arising in the formula for form factors can be written as: 

Qo = i, 
Qi = l, 

Qn = cr 1 cr n - 1 P n , n>3, (1.7) 

where 

P n = det n _ 3 . (1.8) 
Here Ey is an (n — 3) by (n — 3) matrix with the entries 

Ejj = <7 3 i_2j+i, 1 < «, j < n - 3. (1.9) 
The index n — 3 in the expression det n _3 denotes the dimension of the matrix E^. 
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After the Wick rotation to the Euclidean space, a correlation function of the operator 
G can be presented as an infinite series of form factors contributions 

d n f3 



(e(ar)6(0)> 



E 

n=0 



n!(27r) 



■<vac|e(o)|A, . . .,p n )(p n , ■ • • , A|e(0)|vac> 



x exp 



n 



-mr cosh j3j 

3=1 



[1.10) 



where r 

In the present paper we sum up this series explicitly. Now let us present a plan of the 
paper. Section 2 is devoted to a transformation of the determinants to a form, which is 
convenient for summation. In section 3 we introduce auxiliary quantum operators — dual 
fields — in order to factorize an expression for a correlation function and to represent it 
in a form similar to "free fermionic case" . In section 4 we sum up the series ( 1.10 ) into 
a Fredholm determinant. In section 5 we use the Fredholm determinant representation 
for derivation of an asymptotic behavior of the correlation function. 



2 A transformation of the form factor 



A determinant of a linear integral operator I + V can be written as 

f V(xx,xx) ■ ■ ■ V(xx,x n )\ 
r dxx ■ " dx n d . V(x 2 ,xx) ■ ■ ■V(x 2 ,x n ) 

\V(x n ,xx) ■ ■■V(x n ,x n )J 



det(I + V) = 



(2-1) 



Thus, in order to obtain a determinant representation for correlation functions one need 
to represent the form factor expansion ( 1.10|) in the form ( |2.1|) . Determinants of integral 
operators, which we consider also can be called Fredholm determinants. 

The form factor ( |1.3| ) is proportional to the polynomials Q n ( |1 . 7| ) . The polynomial 
Q n for n > 3 is proportional to the determinant of the matrix S^- ( [L.9|) . 

Let M be the following n by n matrix: 



M 



o- 3i _2j_i, 1 < i,j < n. 



Note that M 



j = Sx,j, M 2j = dxjCTs + 5 2 ,jcrx and M nj = 5 n ja n -x for j = 1, 

ij 



(2.2) 
n. Using 



the relation M~ = E,_ 2 i-i (i, 7=3 



2,i-2 {I'iJ — <$, . . . , n — 1), it is possible to show that 



Q, 



det„M, 



n > 0. 



(2.3) 
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The matrix My consists of n 2 different functions, depending on the same set of arguments 
X\ , . . . , x n . 

The main goal of this and next section is to transform the matrix ( |2.2| ) to such a form, 
that entries of a new matrix would be parameterized by a single function, depending on 
different sets of variables, (like V(xi,Xj) in flOD ) 

Mij -> L>ij , D i:j = D(xi,Xj). (2.4) 

In order to study correlation functions, we need to find a square of the polynomials: 

Q£ = det n (Cii), (2-5) 

where 

n n 

Cjk = (M M)j k = ^ °"3i-2i-lC r 3t-2fc-l = X! Cr 3i-2j-lO"3i-2A:-l- (2.6) 
i=l i=— oo 

Here we used the relation = for k < 0. 

Note that the elementary symmetric polynomials have the following expression: 

Z7U J-y Z m=1 

where an integration contour 7 is a circle around origin in positive direction. 

Substituting the above expression into (p.6|) and summing up an infinite series, we 

haVe / dzi r dz 2 z 2n ~ 2j+1 zl n ~ 2k+1 n 

C Jk=f tt 1 f 77^— { vf 1 J{{zi+x m ){z 2 + x m ). (2.8) 

Here we have chosen the radius of the circle 7$ (z = 1, 2) to be greater than one in order 
for the series to converge. 

The matrix C still depends on n 2 different functions Cjk- However, this matrix can 
be transformed to a more convenient form. Let us introduce the following matrix 



1 d n ~^ n 

A]k = ¥^]y- MS x2 + ^ 



(2.9) 

c 2 =0 



with a determinant 

n 

detA = H(x 2 -x 2 ). (2.10) 

i<j 

Using this matrix, we define another matrix D, which differs from C by linear transfor- 
mation 

D = A T CA. (2.11) 



We have an explicit expression for matrix elements of D: 

Djk = j2^iyf d 2 z^^^Y{z 1 ,x j )Y{z 2 ,x k ), (2.12) 



where 

z 2 + x 

and 



Y(z,x) = ^±L i , (2.13) 



J(z) = U(z + x a )(z* + xl). (2.14) 

a=l 

Taking the integral with respect to z 2 , we have after the symmetrization of the integrand 
D jk = — — - - d) ^{^u l Y{u- l z, Xj )){j: u™Y{u-™z-\x k )). (2.15) 

Here u = e 2m ^ and the integration contour 7 is a circle whose radius is larger than one. 
Explicitly, the sum in / can be written as 

3 

cu l Y(cu- l z, x) = Y(z, x) + ujY(uj- 1 z, x) + uj~ 1 Y(ojz, x). (2.16) 

1=1 

Determinants of matrix C and D are related by 

n 

det n C = H(x 2 - x 2 )- 2 det n D. (2.17) 

i<j 

Thus, we obtain a determinant representation for Q\: 

Ql = n detnD • (2.18) 

n - * 2 ) 2 

3 Dual fields 

Entries of the matrix Dj k are parameterized now by a single function D. However, 
an element Dj k , is still not a function of two arguments only, because of the product 
J( z ) — Um=i( z + x m){z 2 + x 2 m ) . This product depends on all x m . In order to get rid of 
these products we introduce auxiliary Fock space and auxiliary quantum operators — dual 
fields. 

Let us define 

$i(x) = q^x) +p 2 (x), ®2{x) = qi{x) +Pi(x), (3.1) 
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where operators pj (x) and qj (x) act in the canonical Bose Fock space in a following way 

(0\ qj (x) = 0, Pj (x)\0) = 0. (3.2) 
Non-zero commutation relations are given by 



\pi(x),qx(y)] = [P2(x),q 2 (y)} =£(x,y) = log((x + y)(x 2 + y 2 



(3.3) 



Due to the symmetry of the function £(x,y) = £(y,x), all fields commute with 

each other 

[* j (x),$ k (y)] = Q, j, k = 1,2. (3.4) 

Dual fields are linear combinations of canonical Bose fields , see page 210 of ||10| . 
Instead of Y(z,x) let us define an operator valued function 

e *i(«) 



Y{z,x) 

Z 2 + X 2 

Instead of Dj k (|2.15 ), we shall introduce an operator in the auxiliary Fock space: 

1 n J y 3 3 

b[x ^ = f -(T,^- l z,x))(Y: uj m Y(uj- m z~\y)). 



(3.5) 



(3.6) 



It is easy to show that an exponent of dual field acts like a shift operator. Namely, if 
g ($i(y)) is a function of $i(y) then 



(oi n 



3 *2 (Xm) 



g ($x(y)) |0) = (0\g q x {y) + £ £(x m ,y) |0) = <?(log J(y)). 



m=l 



m=l 



Using this property of dual fields one can remove the products J(z) from the matrix Dj k . 
For more detailed derivation one should look in at formula (3.6) of the paper ||. 
Standard arguments of quantum field theory show that 



det n -D = (0|det n (^D(xj, X/Je 2 



4 $2 (Xj)+ 3*2(3*;) 



|0). 



(3.7) 



Up to now, we wrote Q 2 n factor of |-F n | 2 as a determinant. An absolute value of the 



form factor is equal to 



l-^nCAb • • • ) AO 



\H n \ 2 Qlll 



i<j 



/(A - A 



\Xi -\- Xj 



|if n | 2 det n L>n 

i<j 



/(A - A 



(X? - £j)(£i + Xj) 



(3i 
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In order to factorize the double product part, we introduce another dual field 



$o0) = %{x) +po{x). 



(3.9) 



As usual 

(0|g (x) = 0, p (x)|0)=0. (3.10) 

Operators qo(x) and po(y) commute with all pj and qj (j = 1,2). The only non-zero 
commutation relation is 

\po(x),q (y)\ = V(x,y), (3.11) 



where 



T]{x,y) = r)(y,x) = 21oj 



(3.12) 



(x 2 -y 2 )(x + y) 

Here we have used the fact that \f({3)\ is a symmetric function : \f(—/3)\ = \f({3)\ ■ It 
is worth mentioning also that the right hand side of ( |3 . 1 2| ) has no singularity at x — y, 
because /(/?) has the first order zero at (3 = 0. Using some of equations for minimal form 
factor from Q 

sinh (3 



/(/3)/(/3 + ztt) 
and f(in) = 4v 2 (0), we can see that 

f(0) 

Hence 



sinh (3 — i sin(7r/3) 



rj(x, x) 



2>/3v 2 (0) ' 
-2 log | Ax 3 1 



where 



A = 8V3w 2 (0). 
Newly introduced dual field also commute 

[%(x),%(y)} = 0=[%(x),<S> J 

Due to the Campbell-Hausdorff formula, we have 



(oi n e 4 °^)|o) = n e^ x ^ = \- n n x- 3 n 

m=l *)i=l m=l i<j 



(3.13) 

(3.14) 

(3.15) 
(3.16) 

(3.17) 



/(log f^) 



[X 2 - X 2 ){Xi +Xj) 



(3.18) 
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Combining these results, we can represent a square of an absolute value of the form factor 
clS cL determinant: 



iWi,...,Aor 




12 n (0|det n (xf\l /2 D(x j ,x k )e^°^ +1 ^ ^) |0). 



Here 



<$>o(x) = $ (x) + $ 2 (x). 



(3.19) 



So we managed to represent a square of an absolute value of the form factor as a deter- 
minant, similar to one of the terms in the right hand side of (|2.1| ). In the next section we 
shall sum up all contributions of the form factors and obtain a determinant representation 
for the correlation function. 

4 The determinant representation 

Because the scaling Lee- Yang model is a non-unitary theory, the normalization constant 
H n ( |1.4| ) is pure imaginary for n odd. This leads to the following relation: 



(vac|6(0)|A, . . .,p n )(J3 n , . . . ,A|e(0)|vac) = {-l) n \F n {(3 u . . . ,/3 n )| 2 . 



(4.1) 



Then the correlation function of G can be written as 



(0(*)0(O)> 

oo „ 

E(-W 

n=0 



n!(2vr) 



d n (3 



n 




xdet n (x.XfcD^^^e^^^^+^^^^e-^ 9 ^)^^))) |0). 



where 




Now we can use (|2.1|) in order to sum up. Finally we obtained a determinant represen- 
tation in terms of an integral operator (Fredholm determinant) 




(4.3) 



S 



The operator U act on a function g(x) as 



[Ug](x) 



dy U(x,y)g(y). 



The kernel of the integral operator U(x, y) is equal to 



U(x,y) = -x 2/J D(x,i/)e5(^W+ ,I> «fe))e-5( e W+ e (f)) 

7T 

where D is given by 

D(x, y) = J-^ f ^ (£ ^(w-'z, *)) ( £ ^(aT™*" 1 , y) 



(4.4) 



(4.5) 



where 



(4.6) 



(4.7) 



z 2 + x 2 

The dual fields $o(x) and $i(x) were defined in the section 3 (see ( [3.1|) and (|3.9| )). The 
main property of these dual fields is, that they commute with each other, so the Fredholm 
determinant det(J — U) is well defined. Certainly det(J — U) is an operator in auxiliary 
Fock space, but it belongs to the Abelian sub algebra. On the other hand, the vacuum 
expectation value of these operators is non-trivial. It follows from commutation relations 



( |3. 3Q , (|3.11| ), that in order to calculate the vacuum expectation value, one should use the 
following prescription 



Mi M 2 Mi Mi Mi M 2 

(oi n e^ a) n e^^io) = n n e^-^ n n e^> x »\ 

a=l b=l a=l 6=1 a=l b=l 



Here 



and 



r)(x, y) = 2 lo£ 



(x 2 — y 2 )(x + y) 



£(x, y) = log((x + y)(x 2 + y 1 



(4.9) 



(4.10) 



Similar Fredholm determinant representations were useful in the theory of correlation 



functions [10, 12, 15, 13 



5 Large r-asymptotic 



Using the technique of [Q, we calculate the long distance asymptotic of the correlation 
function. 
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The determinant can be written as 

det(7 - U(x, y)) = det(I - U(z u z 2 )), (5.1) 



where 



and 



_ poo 

U(z u z 2 ) = dxP 1 (z 1 ,x)P 2 (z 2 ,x), (5.2) 
Jo 

Pi(z,x) = ^- (j2JY(uj- l z,x)j e *<M*H'M (5.3) 

P 2 (z,y) = y (^JY{uj- l z~\y)^ e t^H%) . (5.4) 

The integral operator U[z\,z 2 ) acts on a function g(z) as 

[Uglfa) = iu(z 1 ,z 2 )g(z 2 )dz 2 . (5.5) 



Here the integration contour is a circle around zero in positive direction. 

In the limit r —>■ oo, we evaluate the integral ( |5.2|) by a saddle point method. The 
saddle point of the function 8(x) is x — 1. Hence, we can estimate the integral in ( |5.2j ) 
as 

U(z 1 ,z 2 )=P 1 (z 1 ,l)P 2 (z 2 ,l) f ,& + 0(r- 3 / 2 ) | . (5.6) 

Thus, for the large r asymptotics the kernel U(z\, z 2 ) becomes a one-dimensional projec- 
tor, and its Fredholm determinant is equal to 



det(J-C/)^l- fdzU(z,z). (5.7) 
Using the commutation relations of dual fields, we can evaluate (0\U(z, z)\0) as follows: 



(0\U(z, z)\0) = ——J—Uz^z- 1 ) + ..., (5.8) 
S^tAz V Tar 

where 

3 

Y 1 (z) = u\oj- 1 z + 1) = 3z. (5.9) 
i=i 

The above result leads to the asymptotic form of the vacuum mean value of the Fredholm 
determinant 

/7T — TYIV 

(0|det(J - U)\0) = 1 - \ J + . . . , (5.10) 
which agrees with the large distance contribution from zero- and one-particle states M. 
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Summary 



We considered scaling Lee- Yang Model and we obtained determinant representation for 
the correlation function of the trace of the energy-momentum tensor. We believe that 
determinant representations is the universal language for the description of quantum 
correlation functions in the integrable (exactly solvable) models of quantum field theory. 
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